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I.  INTRODUCTION 

A.   COHERENT  CORRELATION  RECEIVER 

The  optimum  receiver  which  will  detect  the  presence  of  a 
signal  or  discriminate  between  two  different  signals  in  the 
presence  of  additive  white  Gaussian  noise,  is  the  well-known 
coherent  correlator  receiver  structure  which  has  two  equiva- 
lent forms,  as  shown  in  Figure  6.1  [Ref.  2]. 

The  binary  communication  problem  is  modeled  using 
hypothesis  testing  theoretic  principles  in  which  one  of  two 
signals,  S, (t)  or  S0(t),  is  received  in  the  time  interval 
[0,T].  Due  to  the  presence  of  the  noise,  the  observable 
signal  r(t)  takes  on  one  of  the  two  following  forms 


H,:  r(t)  =  S,(t)  +  n(t) 


O^t  ^T 


or 


HB:  r(t)  =  S0(t)  +  n(t) 


0£t  <T 


where  the  noise  n(t)  is  assumed  to  be  a  sample  function  of  a 
stationary  white  Gaussian  process  having  zero  mean  and  power 
spectral  density  (PSD)  level  N0/2  watts/Hz. 

The  optimum  receiver  (in  the  sense  of  minimum  prob- 
ability error,  Pe)  generates  the  test  statistic  G,  where  G 
is  given  by 


(r,Sj)  +  -± 


HSolf 


or  equivalently 


G  =  (r,S,)  -  (r,S0)  +  j- 


and  compares  it  to  a  threshold  in  order  to  decide  whether 
signal  Si(t)  or  S0(t)  has  been  transmitted  in  the  time 
interval  [0,T].   The  inner  product  notation  (•,•)  implies 

(x,y)  =  fx(t)y(t)  dt 

Jo 


and  the  norm  notation 


indicates 


||x||  =  (x,x)  =  fV(t)  dt 


In  order   to  compute  receiver   probability  of  error,    it  is 
necessary  to  determine  the  probability  density  function  of  G 


conditioned  on   the  hypotheses 


H  ,  and 


H 


o  • 


Denoting   these 


density  functions  P(  (G)  and  P0 (G)  depending  on  whether  S,(t) 
or  S0(t)  respectively  was  transmitted,  the  error  probability 
Pe  can  be  expressed  as 

Pe  =  P(H,)|  P,  (G)  dG  +  P(HJ  Jp.  (G)  dG 

loo  '* 

where 

x=  Jk  9  -p(Ho)- 

2-  /71  P(H,) 

and  P(H|)  denotes  the  priori  probability  that  the  hypothesis 
H|  is  true.   It  can  readily  be  shown  that  if  P(H,)  =  1/2, 


_Z 


(1.1) 


e    d 


z 


where 


E* 


-p)e/n0] 


-o-jfsfct)  +  S*(t)J       dt 


and 


Pi 


s.(t)S0(t)  dt 
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We  can  interpret  E  as  the  average  energy  per  bit  and  p  as 
the  normalized  signal  cross  correlation.  Using  the 
following  definition  of  the  error  function 

"°         -l/l 

dz 

Equation    1.1    can  be   written   as 


erf^(x)  "Lme 


F>    =  erfc, 


y<l  -0)E/Uc 


From  the  above  equation  it  can  be  noted  that  as  (l-p)E/N0 
increases,  the  error  probability  decreases.  For  fixed  E/N0 
therefore,  the  optimum  signals  choice  is  one  for  which  P  = 
-1.  This  occurs  whenever  S,(t)  =  -S0(t).  A  system  employing 
this  choice  of  signal  is  known  as  an  optimum  antipodal 
signaling  binary  communication  system. 

B.   JAMMING  OF  COHERENT  RECEIVERS 

1 .   General 

The  objective  of  this  section  is  to  summarize  previously 
derived  results  on  the  performance  of  the  coherent  binary 
(optimum)  receiver  of  Figure  6.1  operating  in  a  jamming 
environment.  That  is,  the  transmitted  signal  is  interfered 
by  the  presence  of  a  jamming  waveform  as  well  as  additive 
white  Gaussian  noise  [Ref.  1] .  The  signal  appearing  at  the 
front-end  of  the  receiver  can  be  mathematically  modeled  by 


H.  :  r(t)  =  S,(t)  +  n(t)  +  n.(t) 

or 

H„  :  r(t)  =  Sa(t)  +  n(t)  +  n.(t) 


CKt  ^T 
0  £t^T 


(1.2) 


where  n.(t)   is  a  jammer   waveform  modeled  as  deterministic, 
yet  unknown  to  the  receiver. 


11 


2 .   Jamm^ r  Op t  im izat ion 

The  effect  of  a  jammer  waveform  on  the  optimum  cohe- 
rent receiver  is  now  presented  in  terms  of  its  impact  on  the 
receiver  probability  of  error. 

In   order   to   determine    the  performance   of   the 

receiver,  we  analyze  the  test  statistic  G,  which  is  given  by 


+  7TII  at     -  Ks" 


with  r(t)  now  given  by  Equation  1.2  .  Due  to  the  assumption 
of  a  Gaussian  zero  mean  noise  having  PSD  level  N0/2 
watts/Hz,  it  is  simple  to  show  that  G  is  a  conditionally 
Gaussian  random  variable  with  mean  value 


E  (G|S.J  =  (S;,Sd)  +  (n:,Sd)  +^[||SJ2 
and  variance 

Var{G|SA}  =-^||Sjf  i  =  0,1 


l|S,|f 


,  i  =  0,1 


With  the  mean  and  variance  known,   the  resulting  conditional 
density  function   of  G  is  given  by 


exp 


g  - 


/27cTS|Sj|* 

V    2  II  J" 

(S^.Sj)  -  (n.,S.)  - 


INI*  -    IIS, |f 


2  1  Vz 


2-i^  .  Is 


k   r  o, 


From  the  general  probability  of  error  expression  (Equation 
1.1)  and  an  approriate  change  of  variables  in  the  integrals, 
when  P(H,)  =  P(H0)  =  1/2  the  expression  for  Pe  becomes 

-~Sa  [Nj  +  d 


B 


/ 

2 


CO 


-tfz 


J1 
I 


'].,    1 


(1.3) 


S^-d] 


-CO 


J 
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where 


Nj   = 


2_ 

%  II  s| 

II  Sail* 


!£. 


2 


d   -    (n.,Sj) 


From  the  jammer  standpoint,  the  optimum  jammer  waveform  must 
be  chosen  in  such  a  way  that  it  maximizes  the  receiver  prob- 
ability of  error.  A  determination  of  the  optimum  jammer 
waveform  can  be  made  by  first  evaluating  the  derivative  of 
Pe  with  respect  to  d,  which  is  the  cross  correlation  between 
the  jammer  waveform  and  the  signal  difference  Sj(t).  Based 
on  Equation  1.3,  we  can  show  that 

L       /si 


8Pe_ 


e 


;inh 


^NJ 


ad  J2TT 
Due  to   the  behavior   of  the   hyperbolic  sine   function,   we 

observe  that 


3d  " 


>0 
=  0 

<0 


,  d  >0 
,  d  =  0 

,  d  <0 


and 


9Pe 


3d 


2  - 


J±Jk 

/2tT 


exp 


-  S^  (d  + 


«$h  ]  > 


o 


for  all  values  of  d.  Since  Pe  is.  a  monotonic  function  of  d, 
it  is  apparent  that  making  d  as  large  as  possible  in  magni- 
tude results  in  the  largest  possible  increase  in  Pe .  In  the 
limit  as  |d|-->co  ,  it  is  easily  seen  that  Pe  -->  1/2. 
However,  from  the  Cauchy  -  Schwarz   inequality 
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(n;>sd>l  <H'IISJ 


(1.4) 


with  equality  if  n.(t)   is  proportional  to  S^ ( t ) 


gy 


I  |n.  |  |  =  /|2     , where   P^.   is  the   jammer  ener 

j     »  '  I;  j 

|  d  |  —  >  oo   implies  that  Pn.   -  -  >  oo        when  |  |  Sj  |  |  <  co 

•J 


Defining 

the  term 

Since 

it  is  not  possible  to  have  infinite  jammer  energy,  one  must 
constrain  n.(t)  such  that  P,,.  is  finite.  From  the 
Cauchy-Schwarz  inequality,  Equation  1.4  can  be  made  into  an 
equality  if 


n.(t)  =  k.Sd(t) 


where  K  is  a  constant  of  proportionality.  If  |  |n.|  |  =  /PM.  , 

I -J        j 

K  must  be  set  to  the  value  7  P^.  / |  1 S j |  )  .  Thus  |d|  is  maxim- 

ized  by  setting 


nJCt)  =/^Tsd(t)/  ||S, 


(1.5) 


and  from  the  above  discussion,  this  results  in  Pe  being 
maximized  also.  Substituting  Equation  1.5  into  the  prob- 
ability of  error  expression  (Equation  1.3)  yields 


Pe  = 


ISrf 


dx 


<i  N.  I   g        yjP      j 


-co 


(1.6) 


It  is  possible  to  put  Equation  1.6  in  a  more  meaningful  form 
by  defining 

E  * 

N   -  SMR   :   signal  to  noise  ratio 
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6, 


*-    -  JSR      :    jammer  to  signal  ratio 


and    observing    that 


IN 


[s,(t)   -  S0(t)]   dt  =   2  E  (1  -p) 


J 


Then,  the  probability  of  error  (Equation  1.6)  becomes 


r 


■CO 


-Jsnr(Ji-P  t  v^jsr) 


R>    =  7 


e 


-  x*/2 


/27C        L.  dx  +//2^T      <2 


-x*/2 


dx 


(1.7) 


Analysis  of  Equation  1.7  indicates  the  fact  that  for  a  value 
of  JSR  beyond  a  so-called  'break  point',  that  is  JSR  > 
(l~P)/2,  an  increasing  SNR  causes  Pe  to  become  worse,  i.e., 
it  increases . 

3.   Optimum  Jammer  Waveforms  for  P SK ,  FSK  and  ASK 

The  effectiveness  of  a  deterministic  jammer  waveform  on 
various  modulation  techniques  will  now  be  presented  for  PSK 
modulation,  where 

S,(t)  =  A  sinWct   ;  S0(t)  =  A  sin(Wct  +  K)        O^t^T 


with  the  constraint  that  WcT  =   n  j\_    \      n  an  integer 
results  for  FSK  modulation  will  be  shown,  where 


Also  , 


S,(t)   =  A  sin  W,t        ;   S0(t)   =  A  sin  W0t 


O^t  ±T 
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with  the  constraint  that   (W,  -  W0)T  =  1 71    \  (W,  +  W0)T  = 

m7L  *,  1  and  m  integers.  Finally  for  ASK  modulation 

S,(t)=A,S(t)   ;   S0(t)  =  A0S(t)    O^t^T 

where  we  assume  that  | |S| |  <  co     and  for  convenience,  that  A( 
>  A0.    The  optimum  jammer  waveforms  and  resulting  receiver 
performances  can   be  obtained  using   Equations  1.5   and  1.7. 
Thus  the  optimum  jammer  for  PSK  is  given  by 

2 


nAt)-JPnl-^r     cos  Wct     0  ^t  £  T 


Since  p  =  -1  for  PSK,  we  have 


P  =  — 


r 


-  xz/2 


dx  + 


J 


-/Sir  (i  +  J~^) 

'TlL       dx 


2SNR  (1  V-JSR) 


-co 


For  FSK  modulation,  the  optimum  jammer  is  given 
nj(t)  =  /ptij  J=   aan^CW,  -WJtcosi  (W,  +  W„  )t 

V  1 
Using  the  assumptions  on  W(  and  W0  ,  the  value  of  /0  is  zero. 
Thus  Equation  1.7  becomes 


1 


O 


-  x2"/? 


-,/§NR 


^«R  (1  +  /2JSR  ) 


7C  dx     +     /y2  7U  fc- 

y^  (i  -  ^jsr  )   -°° 


-x*/2 


dx 


Finally,  for  ASK  modulation,  the  optimum  jammer  is  given  by 

S(t) 
n.(t)  =  /P, 


and  the  probability  of  error  becomes 


16 


s 


f 


where   now 

P- 

and 


dx 


-  Jim  (v^     +  V2JSR) 


TC 


-xa/2 


dx 


S"R  (/^-/2JSR) 


-co 


2A,A0 
At   +  At 


(A,    -  Afl) 


Af  +  Ao 


a--  i  -p  = 

The  'break  point'  for  ASK  occurs  at  JSR  =  0^/2.  Since  ^  < 
1,  in  terms  of  'break  point'  efficiency,  PSK  is  highest  with 
a  'break  point'  occuring  at  a  JSR  of  1.  FSK  is  next  highest 
in  efficiency  with  a  'break  point'  occuring  at  a  JSR  of  1/2, 
and  ASK  is  lowest  in  efficiency  with  a  'break  point' 
occuring  at  JSR  <  1/2. 


C.   INCOHERENT  RECEIVERS 


1.   General 


In  the  coherent  systems  considered,  the  information 
bearing  signals  were  assumed  to  be  known  exactly  at  the 
receiver . 

In  noncoherent  systems,  the  phase  of  the  carrier 
signals  is  not  available  at  the  receiver  so  that  the  phase 
is  treated  as  a  random  variable  uniformly  distributed  over 
[0,27C].  As  such,  we  may  expect  the  performance  of  an 
incoherent  system  to  be  degraded  in  comparison  to  the 
performance  of  the  corresponding  coherent  system.  However, 
because  of  their  simplicity,  incoherent  systems  are  widely 
used  in  many  applications. 
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2.   ASK  (On  -  Off  keying) 

In  the  case  of  incoherent  ASK  (On  -  Off  keying),  the 
signals  at  the  front-end  of  the  receiver  in  the  presence  of 
noise  only  can  be  mathematicaly  modeled  by  either 

r>  /   /  (1-8) 

H,  :  r(t)  =  A  sin  (W<t  +  Q)  +  n(t)    O^t^T 

or 

H0  :  r(t)  =  n(t)  O^t  ^T 


The  amplitude  A,  the  frequency  Wc  and  the  time  of  arrival 
are  assumed  to  be  known  except  that  the  phase  (B  is  modeled 
as  a  random  variable  having  an  a  priori  density  function 


0^6  £27^ 

otherwise 


27c 

%>(8)  =1 0. 


The  additive  noise  is  again  assumed  to  be  zero  mean,  white 
Gaussian,  with  power  spectral  density  level  N0 / 2  watts/Hz. 

In  the  absence  of  a  jamming  waveform  n.(t),  the 
optimum  receiver  for  decoding  the  binary  information  (which 
is  transmitted  via  the  use  of  either  a  sinusoid  that  is 
incoherently  received,  or  no  signal  at  all)  in  each  interval 
of  duration  T  sec,  is  the  well-known  quadrature  receiver. 
Its  structure  is  shown  in  Figure  6.2  [Ref.  2]. 

There  is  an  alternate  form  of  the  quadrature 
receiver  obtained  by  replacing  the  correlators  with  matched 
filters  having  an  impulse  response  given  by  h(t)  =  sin  Wc(T 
-  t)  and  h(t)  =  cos  Wc(T  -  t),  0  <  t  <  T,  and  sampling  the 
outputs  at  t  =  T  (Figure  6.3)  [Ref.  2]. 

Still  another  important  alternate  form  of  the  quad- 
rature receiver  is  an  incoherent   matched  filter  followed  by 
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an  envelope   detector  and  a  sampler   as  shown  in   Figure  6.4 
[Ref.  2]. 

3.   FSK 

For   binary   frequency   shift    keying   (FSK),  the 

received  signals   at  the   front  end  of   the  receiver   in  the 

presence   of  noise   only  can   be   mathematically  modeled  by 
either 


H,  :  r 


(t)  =  A  sin  (W,t  +0)  +  n(t)    O^t^T 


or 


H0  :  r(t)  =  A  sin  (W0t  +  Cj)  )  +  n(t)   ■  0  ^t  ^  T 


where  the  additive  noise  is  again  assumed  to  be  zero  mean, 
white  Gaussian,  with  power  spectral  density  level  N0/2 
watts/Hz.  Usually  the  frequencies  W(  and  W0  differ  substan- 
tially. The  phases  0  and  (h  are  statistically  independent 
random  variables,    uniformly  distributed  over   the  interval 

[0,27H]. 

Insofar  as  the  signal  amplitudes  is  concerned,  two 
cases  are  distinguished.  In  the  first  case,  the  information 
bearing  signal  amplitudes  are  known  and  equal.  In  the 
second  case,  the  information  bearing  signal  amplitudes  are 
treated  as  independent  random  variables  that  are  Rayleigh 
distributed  due  to  multipath  propagation.  For  this  second 
case,  the  following  hypothesis  testing  theoretic  model  is 
used 

H,  :  r(t)  =  A  sin  (W,t  +  0)  +  n(t)    0  ^t  ^T 
and 

H0  :  r(t)  =  B  sin  (W0t  +  0  )  +  n(t)    O^t^T 
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where  we  assume  that  the  signal  amplitudes  even  though 
random,  remain  constant  over  a  . T  -  second  interval.  The 
noise  is  modeled  as  in  the  previous  analysis,  and  the  signal 
amplitudes  have  (a  priori)  probability  density  functions 
given  by 


2. 

a 


fA(a)  =  -J5j  exp  (-  75)       3.  >  0 

and 

b         ,* 
fB(b)  =^iexp(-  -fjz)  b>0 

where  A2,   =  E{A2 }  =  E{B2} . 

The  results  associated  with  ASK  incoherent 
receivers  can  be  applied  toward  obtaining  the  optimum 
incoherent  FSK  receiver.  This  is  worked  out  in  Reference  2. 
The  resulting  receiver  structure  is  shown  in  Figure  6.5  . 
Individual  channel  matched  filters  shown  in  Figure  6.5  can 
be  essentially  replaced  with  tuned  bandpass  filters  having 
center  frequencies  corresponding  to  the  "mark"  and  "space" 
frequencies  Wt  and  Wo. 
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II.  ANALYSIS  OF  JAMMING  ON  ASK 

A.   GENERAL 

In  this  chapter,  the  effect  of  a  deterministic  jammer 
waveform  on  the  performance  of  an  incoherent  ASK  receiver 
will  be  investigated. 

The  signal  at  the  front-end  of  the  receiver  is  given  by 
Equation  1.8  with  the  modification  that  under  either 
hypothesis,  a  jamming  waveform  n;(t)  is  present  during  the 
time  interval  [0,T]. 

For  the  coherent  ASK  receiver,  the  optimum  energy 
constrained  jammer  waveform  could  be  obtained  by  using 
Equation  1.5  .  However,  the  complexity  of  the  mathematical 
expression  for  probability  of  error  in  the  noncoherent  case 
makes  it  very  difficult  if  not  impossible  to  derive  the 
optimum  energy  constrained  jammer  waveform  in  closed  form  as 
will  be  seen  in  the  sequel. 

A  reasonable  postulation  is  that  the  optimum  jamming 
waveform  for  the  coherent  ASK  receiver  can  act  as  a  good 
jammer  for  the  incoherent  ASK  receiver  also.  Thus,  such 
near  optimum  jammer  signals  are  studied  and  evaluated  in 
terms  of  their  effect  on  the  performance  of  the  noncoherent 
ASK  receiver. 

B.   ANALYSIS  WITH  NEAR  OPTIMUM  JAMMER 

Analysis  of  the  incoherent  receiver  starts  from  the 
mathematical  model  of  the  receiver  front-end  input  waveform 
r(t)  given  by  either 

H,  :  r(t)  =  A  sin  (Wct  +  Q  )  +  n(t)  +  n.(t)    0  £t  IT 


■j 


or 
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H0  :  r(t)  =  n(t)  +  n.(t)  O^t^T 

j 

where  n.(t)  is  the  jammer  waveform  present  during  the  time 
interval  [0,T] . 

In  the  absence  of  n.(t),  the  optimum  receiver  for  the 
binary  ASK  problem  is  well-known.  Its  derivation  is  well 
documented  in  the  statistical  detection  theory  literature 
[Ref.  2].   The  receiver  structure  is  shown  in  Figure  6.2 

In  this  section,  the  effect  of  n.(t)  on  this  receiver  is 
analyzed  by  evaluating  the  resulting  Pe  ,  under  the  assump- 
tion that  n.(t)  is  a  deterministic  jammer  waveform,   however 

sJ 

unknown  to  the  receiver  itself.  Receiver  performance 
requires  determination  of  the  statistics  of  either  G2  or  G, 
where  G2  is  the  output  of  the  quadrature  detector  and  is 
given  by 

1  X  3- 

G  =  X  +  Y 

where 


X  =  Jcr(t)  sin  Wct  dt  Z      (r,S) 
Y  =  (Tr(t)  cos  Wct  dt  =   (r,C) 

Provided  that  the  random  variable   £p   is  fixed  to  some  value 
9 ,  X  and  Y  are  conditional  Gaussian  random  variables  with 


E{X|H,  ,8]  =  (ASa,S)  +  (n.,S) 


A 


m>|9 


and 


{y|h,,0}=  (ASe,C)  +  (n;,C)  I  my(e 
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where  Se   is  used  in  place  of   sin(Wc  t  +  Q   ).    It  can  also 
be  seen  that 


Var  {x|H(,QJ  =  E  {(n,S)a} 


N.T 


sin  2WCT 


2Wt_T 

=  Var  {x|HaJ 
and  similarly 


sin  2ty.T 


~l 


>WCT 


=  Var 


{y|h0} 


For  convenience  we  assume  WcT  =  n 7C ,  where  n  is  an  integer. 
Thus  the  sine  terms  above  vanish  resulting  in  var{X  |  |-|.  ,  0  } 
=  var{Y  |  t-{,  ,  ^  }  =  N0T/4  =  (J  2  .  If  this  assumption  is  not 
made,  an  additional  term  results.  However,  if  7C /Wc  <<  T, 
the  additional  sine  term  is  small  and  may  be  neglected. 
Furthermore  the  covariance 


X  -  E  {> 

'IH.,0}1 

• 

Y 

E   |(n,S)    (n,C)j 

N0T 

1  -  cos  2WCT 

4 

_ 

2WtT 

{y|h.,0 


H;.0 


=  0 


i  =  0,1 


provided  the  assumptions  on  Wc  hold.  This  implies  that  for 
any  given  value  of  0  ,  both  X  and  Y  are  uncorrelated 
Gaussian  random  variables  and  therefore  statistically  inde- 
pendent. The  density  function  of  G2  conditioned  on  the 
phase  0  is  noncentral  Chi-Squared  distributed,  and  is  given 
by 
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^c«|H„e,.^.^-ii£!^E) 


,  g>0 


and  zero   otherwise.   The  variance  (J  2    is  defined   above  and 
becomes 


(2.1) 


S  = 


m*ie 


+  m*„  =  EZ   {x|H,,9}  +  E2  {y|H,,g} 
[(AS»,S)  +  (nj,S)f+   f(ASe,C)  +  (n^C)"]* 


Due    to    our   assumption   WcT    =    r\7t  ,    we    obtain 


f  AT  \  r  -i 

S  =  \^2~)  +  AT   [(n.,S)  cos  9  +   (n.  ,C)   sin  0 

+   (n    ,sf    +   (n    ,cf 


Similarly,  if  r(t)  consists  of  noise  and  jammer  only,  then  X 
and  Y  are  also  independent  Gaussian  random  variables  with 


4 

=  m. 


E  {  X|H0|  =  (n.,S) 
and 

e|y|H0|  =  (n.,C)   =  mY 
The  density  function  of  G2  assuming  no  signal  is  sent  is 
V(,|H.,.  _L  ^{.^)lM ) 


and  zero  otherwise,  where 

S'~<     +mr  =  Cn.,Sf  +  (n.,cf 


(2.2) 
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The  resulting  average  probability  of  error  is  given  by 


P*   = 


pp-P(HJ+  PM   P(H.)  x 

P(Ho)  /  V(8lHo  )  dg  +  P(H.)    |fy  (g|H|   )  di 


-00 


where  Tj    is  the  threshold  with  which  G  is  compared  in  order 
to  decide  which  is  the  true  hypothesis.  Observe  that 

fjzCglH.)  =   f^|e  (g|H,,  0)  f9(0)  dQ 

The  first   integral  in   the  above  expression    for  Pe   is  Br 
(probability   of  false   alarm) .     It   can  be   expressed   as 
follows 


■oo 


2o 


x exp  C- 


g  +  s 


n 


2  o2 


-)  I0  (  ^=-  )  d 


g 


V  exp  (-  liil)  i  (  a'v)  dV 


/7?/(J 

Q  (  a',  n/o  ) 


where 


.CD 


Q  (  a  ,  6  )  =     V  exp  (-   V  t  Q 


# 


2 


is  a  well-known  tabulated  function  called  the  Marcum 
Q- function.  The  second  integral  in  the  above  expression  for 
Pe  is  PM  (probability  of  miss).  It  can  be  expressed  as 
follows 


V 


tv-~  r* 


G*|Q(g|  H,  ,8)  fg(8)  d9 


dg 


(2.3) 
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rco 


r   r^J- 


f  (8) 


-co 


J, 


^f   (g|H, ,6)  dg 
G  |  8 


L-  *  0 


de 


where  in  the  second  equality  the  order  of  the  integration 
has  been  changed.  The  inner  integral  of  Equation  2.3  can  be 
expressed  in  term  of  the  Marcum  Q- function  as 


"^i 


2  a; 


exp  (-  £_l_i)  I0(€L)    dg 


2  a 


(% 


exp  (-   V  +  a  )  lfl(  aV)  dV 


=  1  -  Q  (  a  ,  n  /  a  ) 


where  q(   =  JS//T  .   Therefore  the  probability  of 


error  can  be 


written  as 


/s7 


Pe  =  P(H0)  Q  (— ',  n/o  ) 


(2.4) 


*  +  P(H.) 


f2  7T 


1  - 


2  TT 


y^ 


Q  (~  ,  n  /o  )  di 


so 


where  the  dependence  on  Q  is  imbedded  in  the  term  /  S  and 
the  threshold  Tj  for  the  jammer  absent  case  is  obtained  as 
the  solution  to  the  equation 


-  A  T/2N0 
Q  I0(,2An/N0)  =  R 


PCH0) 

5    R  ~  P(H,  ) 


which  can  be  equivalently  put  in  the  form 


e 


-  A  T/2M0 


h 


A  T/N- (  n  /  a   )    =  R 


The  I0  ( •  )    function  used  here   and  also  previously   used  in 
conjunction   with  the   development   of   Pe  is   the   modified 
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Bessel  function  of  the  first  kind.  The  appropriate  setting 
of  the  threshold  is  an  important  issue  for  both  coherent  as 
well  as  noncoherent  receivers.  One  could  use  the  threshold 
setting  that  would  be  derived  from  the  analysis  of  receivers 
operating  in  additive  white  Gaussian  noise  interference 
only.  This  approach  can  be  quite  unsatisfactory  as  demon- 
strated in  Reference  1.  A  better  approach  would  be  to 
obtain  Pe  as  a  function  of  the  threshold,  and  then  search 
for  the  threshold  that  minimizes  Pe.  While  this  approach  is 
intuitively  appealing,  it  is  often  mathematically  intrac- 
table. While  the  threshold  issue  is  not  addressed  in  this 
particular  section,  it  will  be  discussed  in  more  detail  in 
the  sequel  and  simulation  results  will  be  presented. 

If  the  definition  of  average  signal  energy  previously 
introduced  is  used,  we  have  E  =  A2T/4-  which  is  reduced  in 
half  in  comparison  to  the  binary  signal  transmission  case, 
due  to  the  fact  that  the  information  bearing  signal(s)  do 
not  have  equal  energy.  In  order  to  afford  comparisons  with 
the  coherent  receiver  case,  we  will  implicitely  boost  the 
value  of  signal  amplitude  A,  to  obtain  E  =  A2T/2  in  order  to 
have  agreement  with  previous  cases  insofar  as  signal  ener- 
gies is  concerned.  Thus  the  threshold  determination  equation 
now  becomes 


-SNR 


e 


[\/2SNR  (  n  /a  )] 


=  R 


R  = 


P(H„) 

F(H,  ) 


(2.5) 


If  we  assume  P(H.)  =  P(H.)  =  1/2,  then 


L       -   »J0  u 

and  the  threshold  setting  equation  becomes 

-SNR   f  

Q  I0[/2SNF  (1  /„  )J   =  l 


iri  /a  ) 
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If  we  are  to  find  the  optimum  jammer  waveform  so  as  to 
maximize  Pe,  an  attempt  must  be  made  to  solve 

—  =  0  and  as,  =  0 

Unfortunately  the  resultant  equations  are  mathematically 
involved  and  do  not  appear  readily  solvable  for  n.(t).  It 
seems  however  that  a  'good'  jammer  waveform  can  be  postu- 
lated based  on  the  results  obtained  for  coherent  ASK.  It  was 

found  for  that  case  that  the  optimum  n.(t)   is  a  tone  at  the 

o 
carrier  frequency.  Thus  the  following  jammer  waveform  can  be 

used  as  a  potential  near  optimum  jammer,  namely 


.  /  p  f  2  6  ") 


Observe  that  with   this  choice,   ||n.||2  =   P*.  .    The  prob- 

J  J 

ability  of  error  Pe  can  now  be  determined  using  the 
threshold  setting  equation  (Eqn.  2.5)  and  the  previously 
derived  expressions  for  S  and  S'.  The  effect  of  the  near 
optimum  jammer  waveform  on  the  receiver  (i.e.,  incoherent 
receiver  performance)  can  be  analyzed  by  evaluating  Pe  as  a 
function  of  n.(t)  using  Equation  2.6  .  Note  that  in 
Equation  2.6  the  jammer  energy  is  P^.  .   It  can  be  shown  that 


(n.,C)  =fr/n~/X  sin  Wrt  cos  VLt  dt  =  0 


Then  from  Equations  2.1  and  2.2,  the  probability  of  error  Pe 
given  by  Equation  2.4  can  be  expressed  in  terms  of  SNR 
(signal-to-noise  ratio)  and  JSR  (j ammer- to- signal  ratio) 
using  the  fact  that 
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S  o  / 

~qJ=   2   SNR  (1  +  2  /JSR  cos9+  JSR) 
"tS;  =   2   SNR- JSR 

and    the    fact    that    Tj  /^-         can      be    obtained    from   Equation   2.5. 
For      equally      likely   hypotheses,         the      probability      of 
receiver   error    (Equation   2.4)    can   be    expressed   by 


i-Fl  +J3  (  /2SNRJSR   ,  ^/(J  ) 

__L.(^(    /2SMR  (1  +  2v/JSl?cos9+  JSR   )    ,"7/<j)   d 9J 


(2.7) 


From  Equation  2.7,  Pe  greater  than  0.5  can  occur  if 


Q  (a  ,3  )>^-  Q  (a',6  )  di 


(2.8) 


where 


Ct'=    /2SNR  (1  +  2/JSR  cosO+  JSR) 
#  =  ^2 SNR  JSR    and  ^  =  7?  /<j 


Since  the  periodic  sinusoidal  function  imbedded  in  (A  '  is 
integrated  for  one  period,  the  condition  o(  -  oC  '  might 
cause  the  condition  of  Equation  2.8  to  be  satisfied  for  a 
fixed  value  of  &  .  Thus  there  is  a  possibility  of  obtaining 
Pe  greater  than  0.5  for  a  JSR  value  beyond  the  critical 
value  of  0.25  which  can  be  obtained  from  the  condition 
pi  =  qC  ' ,  that  is, 


J  2 SNR-  JSR  =  y/2 SNR  (1+2  JJCiR   cos0  +  JSR) 


(2.9) 


/JSR  cosB  =  -  — 
2 
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Receiver  performance  in  the  presence  of  a  jammer  is 
expressed  in  terms  of  signal-to-noise  ratio  (SNR)  and 
jammer- to-siganl  ratio  (JSR).  A  rectangular  pulse  of  dura- 
tion T  seconds  has  amplitude  spectrum  AT  sine  Tf,  and  B  = 
2/T  is  a  rough  measure  of  its  bandwidth.  Thus  the  expression 
of  the  form 


A2T/2    AZ/2     Az/z 


(2.10) 


Uo_  2         PSD-B 
2    T 


can  be   interpreted  as   the  ratio  of   signal  power  to  noise 

power  (SNR)  in  the  signal  bandwidth  [Ref.  5].    The  term  JSR 
can  be  expressed  by 

P^.    jammer  energy  (2.11) 


A^T/^  signal  energy 


where  P-n.   represents  the  jammer  energy  defined  before. 

C.   VARIABLE  THRESHOLDING  EFFECT 

It  is  apparent  that  when  a  fixed  threshold  value  is  used 
by  the  receiver,  the  effect  of  the  jammer  waveform  having 
energy  such  that  the  resulting  JSR  value  is  above  a  certain 
level,  is  such  that  the  receiver  may  be  rendered  inoperable. 

It  could  therefore  be  suggested  that  setting  the 
threshold  based  on  the  value  P(H|)  or  P(H0)  only  (see  Eqn. 
2.5),  may  not  be  desireable.  Values  of  threshold  other  than 
some  fixed  value  may  result  in  improved  receiver  perform- 
ance, in  other  words,  reduce  the  jamming  effect. 

Recall  that  the  threshold  setting  equation  is  of  the 
form 

-SNR 


Q  I0(  J2SUR   a)  =  R  (2.12) 
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where  0\  =  ,  l(J  •  we  now  can  attempt  to  reduce  Pe  by  an 
appropriate  choice  of  threshold  J  .  By  varying  the  value  of 
R  (instead  of  using  R  =  1)  we  can  obtain  this  threshold 
which  denote  q(   '  from  the  expression 


72SIIR 


SNR 

i;'(r  e    ) 


where  I0"1  (X)   is  the  inverse  modified  Bessel  function  of  the 

first  kind.    Using  an  approximation  to  I0 (X)  which  is  given 

by 

x 

for  x%>  1  , 


I0(x) 


e 


2  7t  x 
Equation  2.12  can  be  rewritten  as 


(2.13) 


4  yj  2SIIR  a'  -  2  j^a' 

=  2  j2f»2TT   +  J?rK2SNR)  +  USMR  +  fa  R 


From  Equation  2 . 13  ,  it  can  be  recognized  that  when  SNR  is 
large,  variation  of  the  value  of  R  does  not  significantly 
affect  the  threshold  value  qC  '  because  taking  the  logarithm 
of  R  reduces  the  effect  of  R  further.  Thus  the  large  value 
of  SNR  suppresses  the  effect  of  variation  of  the  term  J2n  R 
on  ^'  .  This  limited  variable  thresholding  effect  on  the 
receiver  which  is  under  significant  jamming  is  analyzed  in 
graphical  form  in  Chapter  5  for  various  values  of  R. 
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III.  ANALYSIS  OF  JAMMING  ON  FSK 

A.   GENERAL 

For  binary   incoherent  FSK  with   a  jammer   present,   the 
received  signals  under  the  two  hypotheses  are  either 

H,   :  r(t)  =  A  sin  (W,t  +  0  )  +  n(t)  +  n.(t)     0  *t  ^T 

j  ~ 

or 

W0     :  r(t)  =  A  sin  (W0t  +  (f))  +  n(t)  +  n.(t)     0  It  iT 

By  separating  the  frequencies  W,  and  W0  sufficiently,  we  can 
form  signals  that  are  orthogonal,  have  equal  energy,  and 
have  the  same  advantage  of  ease  of  generation. 

The  modified  FSK  receiver  structure  which  is  capable  of 
varying  the  output  of  the  each  envelope  detector  (which  is 
followed  by  a  multiplier)  is  diagramed  in  Figure  6.6  .  The 
optimum  receiver  for  the  case  where  no  jammer  is  present  can 
be  derived  from  statistical  decision  theory  and  is  shown  in 
Figure  6.5  .  This  receiver  can  be  obtained  by  setting  c(  = 
1/2  in  the  modified  FSK  receiver  shown  in  Figure  6.6  . 

In  practice,  incoherent  FSK  is  widely  used  because  of 
its  simple  receiver  structure,  its  small  performance  penalty 
due  to  lack  of  phase  coherence  and  its  more  efficient  use  of 
signal  energies  in  comparison  to  incoherent  ASK.  In  addi- 
tion, we  are  not  faced  with  the  difficulty  imposed  by  a 
threshold  that  must  change  with  SNR  as  is  the  case  with 
incoherent  ASK.  Thus,  a  receiver  which  is  known  to  be 
optimum  for  incoherent  FSK  transmission,  has  been  modified 
by  including  some  channel  weighting.  This  has  been  done  in 
order  to  be  able  to  determine  whether  or  not  such  channel 
weighting  can  reduce  the  effect  of  the  jammer. 
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This  chapter  is  devoted  to  investigating  the  performance 
of  the  modified  incoherent  FSK  receiver  in  the  presence  of 
jamming  and  additive  white  Gaussian  noise.  By  letting  pC  = 
1/2  we  can  as  a  byproduct  obtain  the  performance  of  the 
conventional  incoherent  FSK  receiver  of  Figure  6.5  in  the 
presence  of  jamming  and  additive  white  Gaussian  noise. 

B.   ANALYSIS  WITH  NEAR  OPTIMUM  JAMMER 

The  modified  receiver  performance  can  be  obtained  by 
introducing  a  null  hypothesis  (no  signal)  as  a  dummy 
hypothesis  and  by  following  the  same  reasoning  as  in  the 
analysis  of  incoherent  ASK  presented  in  the  previous 
chapter. 

The  receiver  function  is  to  compare  the  envelopes  at  the 
output  of  each  channel  once  every  T  seconds  and  decide  in 
favor  of  the  larger  of  the  two  envelopes  (Figure  6.6).  For 
the  purpose  of  analysis,  let  us  first  assume  that  a  'mark' 
signal  has  been  transmitted,  that  is,  the  hypothesis  H,  is 
assumed  to  be  true.  An  error  is  committed  if  V0  exceeds  V| . 
An  error  is  also  committed  if  V|  is  larger  than  Wq  when  a 
'space'  signal  has  been  transmitted,  that  is,  the  hypothesis 
H0  is  assumed  to  be  true  [Ref.  3]. 

Let   Pe,  denote   the  probability   of  the   first  type   of 
error  described  above,  which  is  expressed  as  Prob.(V0  >  V(  | 
H.).   Under   the  assumption   that  a   'mark'  signal   has  been 
sent,   the   output  qf   of  one  of   the  envelope   detectors  is 
given  by 

where 

X,  =   r(t)  sin  W,t  dt   ;   Y,  =   /  r(t)  cos  W,t  dt 
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Observe  that  X(   and  Y,  conditioned  on  the   phase  and  either 
of  the  two  hypotheses  are  Gaussian  random  variables  with 

E  |x,|H,,0}=  /TA  sin  (W,t  +  0  )  sin  W,t  dt 

+  fTn.(t)  sin  W,t  dt  =  (S,,S),  +  (n.,S), 


and 


E  {'Y,|  H,,0}=  /A  sin  (W,t  +  Q  )  cos  W,t  dt 

+  I  n.(t)  cos  W,t  dt  =  (S,  ,C),  +  (n.,C), 


where  S|  represents  the  function  A  sin  (W(  t  +  6  ),  S  repre- 
sents the  function  sin  W| t  and  C  represents  the  function  cos 
W,t.  Likewise,  assuming  again  that  W(T  =  n7L,  n  an  integer, 
and  that  n(t)  is  zero  mean  white  Gaussian  noise  with  PSD 
level  N^/2  watts/Hz,  we  obtain 

Var  [x,|H,,eJ  =  Var  [yJh^G}  =  ^ 
Furthermore,  it  can  be  shown  that 

e  | [x,  -  e  {x,|h,  ,e}][y,  -  eJy.Ih.^Ih^oJ      =  0 

so  that  the  conditional  r.v.'s  X,  and  Y|  are  uncorrelated 
and  therefore  independent.  The  sum  involving  random  vari- 
ables X,  and  Y,,  and  producing  q*,  will  result  in  a  non- 
central  Chi-Squared  distribution  so  that 


f0>fq,|H,,6)  =  ^lexp(- 


^)I„<^> 


2  o 


where 


S„  =  E^XjH,,©]  +  E*{y1|H,,0] 

(S,  ,S),  +  (n.,S),[  +   [(S,  ,C),  +  Cn.,C),] 
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q,  >o 


and  (j  2    -      N0T/4.  Using  standard  random  variable  transforma- 
tion techniques,  it  is  not  difficult  to  show  that 

Wi,iH.,e>  =^«xp(-lL!!)io  fij^L)      q>2i0 


so  that 


(3.1) 


fO,  (cl.  I  H-   >  =    W  (qjH.,6)  %  O)  d0 

-/a  o 


where  the  dependence  on  6  is  imbedded  in  the  term  SM  .  We 
now  need  to  obtain  the  statistics  of  the  output  of  the 
multiplier  following  the  upper  envelope  detector.  (See 
Figure  6.6).  That  is,  the  probability  density  function 
fv  (^i  I  **  I  )  needs  to  be  derived.  From  standard  transforma- 
tion theory,  using  the  relation  V,  =  2(1  -  (/.  )q(,  it  can  be 
shown  that 

f  (v,|H,)  = fQl(  — |H,  ) 

2(1- a)      2(1- a  )  (3.2) 

47t(1.a )  jhill^-  -} 


°    a 


1 2  a 


'{iii^V5'  )de 

a-2- 
On   the  other   hand,   the   output  of   the  lower   envelope 
detector  when  H  |   is  assumed  to  be  the   true  hypothesis,   is 
given  by 
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q„  = 


X  +  Y 


where 


X, 


r(t)  sin  W0t  dt   ; 


Y-  =    r(t)  cos  W„t  dt 


Following  a   similiar  procedure   as  used   above,   it   can  be 
shown  that 


|y0|H(,  e|  =   /A  sin  (W,  t  +  6)  cos  WDt  dt 

+  I  n  .(t)  cos0W  t  dt  =  (S  ,  ,C)0  +  (n.  ,C)0 
Jo 

e|Xo)H,  ,ej  =    A  sin  (W,  t  +  6)  sin  W0t  dt 
Jo 

+      nj(t)  sin0W  t  dt  =  (S,  ,S)0  +  (n.  ,S), 


and 


Here,  S  and  C  have  the  same  meaning  previously  defined 
except  that  the  subscript  "0"  outside  the  inner  products 
implies  that  we  should  interpret  S  as  sin  W0t  and  C  as  cos 
Wct .  It  can  also  be  demonstrated  that 


Var 
and  also  that 


(  Xe|H,,e}=  Var  |yo|H,,G} 


N0T 


e{  [  X0  -  E  [xjH.,0}]  [  Y0  -  E  {Y0|H„eJ]|H,,e}=  0 

so  that  the   conditional  r.v.'s  XD   and   Y0  are  uncorrelated , 
hence   independent.   Thus   similiar   to   Equation  3.1,    the 
expression   for   the   conditional   density   function   of   q 
becomes 
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fQ(qJH. > 


azexp  (- 


qt+  s< 


■2 )  ' 


I  (  <WS°» 


)  d8 


q0^o 


where 


s0,  =  e*{  x0|h,,gJ  +  e*{  y0|h,,q] 

(S,  ,S)0  +  (n.,S)J  +   f(S,  ,C)0  +  (n  C)J 

•J  J  L  O      J 


and   dependence  on   6   now  is   imbedded   in   the  term   S  0 
Applying   again   the   random  variable   transformation  VD 
2(7(q0,  the  conditional  density  function  of  Vo  is  given  by 


fVo(V0|H.)  = 


exp  - 


(^)*  s 


0l 


2a' 


(3.3) 


Vo  __ 
{   2  aySo'  ■ 


-\de 


In  order  to  compute  the  probability  of  error,  we  can  now 
use  the  previous  expressions  for  the  conditional  probability 
density  functions  which  are  derived  assuming  a  'mark'  signal 
has  been  sent.  That  is,  for  a  given  value  of  Vh  an  error 
is  made  if  V0  >  V i .  Thus  the  average  error  probability  is 
found  by  averaging  the  conditional  error  probability  given 
by 


P(V.)  =    f  (V0|H,,V,)  dV0 

A/. 


over    all   Vi  .    That    is  , 


re, 


Prob, 


v0>v, 


I  H,] 


(3.4) 
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-  JTf 

Jo  L  v, 


f„  (V,  H,  )  dV, 


fu  (V,|H,  )  dV, 


V, 


Substituting  for  the  density  functions  in  Equation  3.4,  from 
Equations  3.2  and  3.3,  we  obtain 


roo 


2  7ZJ0    [_J0 


271,     r<° 


oaexp  (  -  ^2a  /  +  b«  \  ' 


(3.5) 


I0  (_2aVf^)  dV0|  de]  fV((V,|H()  dV, 


where  the  order  of  the  integration  has  been  changed  for 
computational  convenience.  In  the  above  equation  the  inner- 
most integral  can  be  expressed  in  terms  of  the  Marcum 
Q-function  as  follows 

—  exp  (-  ^~;  *  So'  )  I.(  2a  Vbo')^Va 


V, 


2a 


=  q  <  ^K 


V, 


) 


a  '  /aa 
Then,  Equation  3.5  for  Pe   becomes 

vS«i 
o 


1 


ei 


ZTT        -^0 


V. 


L     -'O 


rt    '    2aa 


fV((V,|H,)    dV, 


(3.6) 


3*q    (    ^Z  ,     il£L  x   )    di 


J0   L  2  IT     / 


^7C 


I    x    exp    ( - 

L  2tt   J0 


«*♦    Of},    ! .(    x    S 


de 


dx 
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where  in   the  second   equality  the   change  of   variable  x 


V 


,  has  been  made 


2(l-a)a 

From  the  orthogonality  property  of  the  signal  pair  used 
(which  is  obtained  by  assuming  sufficient  separation  between 
two  frequencies  and  that  Wl  as  well  as  W0  are  large),  we 
have 


f 


sin  (W.t  +  9 )  sin  W„t  dt  =  0 


so   that  the   term   S 0|  is   independent   of  Q   .     Therefore 
Equation  3.6  can  be  rewritten  in  the  following  form 


P  = 
e,  2 


Ja      L  Jo 


t  X  +   I  — —  )  \ 

x  exp  (-  2 >  • 


I.  (x 


v  3,i 


)  dx 


d9 


where  the  order  of  integration  has  been  changed. 
Furthermore,  using  the  following  formula  involving  an  inte- 
gral of  a  Marcum  Q-function  [Ref.  A] 

Jn 


<£**? 


1  -  Q  (Jn— 


Vo*+  a; 


o? 


2. 


*   grp^'  Va,^-o-  'Vo^?i 


the  inner  integral  in  Equation  3.7  can  be  simplified  in  such 
a  way  that  Pe,  becomes 

(3.8) 


p  = 


(1-   a) 

(ITaf  +a: 


i    _ 


2tt 


0Q  u     V(l-   af+aa,    a     \/<l-a    f+aa 


d8 
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+   CL, 

(1-a  )+a* 


1 

2~  . 


Q  ( 

0 


a 


a  v(l-a  )+a  ,  a  v(i_a  7+a 


where  the  dependence  on  6     is  imbedded  in  the  term  S  u  . 

Following  exactly  the  same  procedure  used  in  obtaining 
the  expression  for  Pe, ,  it  can  be  established  that  the 
expression  for  Pe0,  which  denotes  the  error  probability  when 
H  0   is  assumed,  that  is,  Prob.(  V!  >  V0|  H0)  takes  the  form 


(3.9) 


a 


(1-  oO+aa 


1  - 


1 
7¥J 


'±n 


Q  < 


(1-  a)   v/5,,  /  O-oO' 


1-a 7+  a1  j    (f  \/   (,_^)>+i<^ 


(1-  a) 

(1-  af+a* 


n-l  o  ( 


/   (1-a  L  \/s^  /    aa    v,a 
V(l-a  7+a*-  ,"5-/(1-  /)V)d? 


where  only  the  term  S  ,»  is   dependent  on  <b  •    Therefore  the 
total  average  probability  of  error  (Pe)  can  be  obtained  from 

Pe  =  -*>   (  Pe,  -  P€e  ) 
assuming  that  the  two  hypotheses  are  equally  likely. 

For  the   special  case  q(    -      1/2,   the  performance   of  the 
optimum  receiver  in  the  presence  of  jamming  is  given  by 


v4»-i 


2tt 


/T"    a    ,/7     a 


(3.10) 


JL  v^ 


(^'T,vf'^\  d8 


2tt 


^ 


f  ^     a     ,  /2~  ~} 

q  (    J.  </§T7    ivg7 
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Observe  that  Equation  3.10  yields  the  performance  of  a 
conventional  incoherent  FSK  receiver  (with  no  channel 
weighting)  in  the  presence  of  jamming.  If  we  now  use  as  the 
jamming  waveform  the  jammer  which  is  optimum  against  a  cohe- 
rent FSK  receiver,  that  is 


I 


nJ(t)  :/p^sinI(W,-WJtcosi(W1  +W0)t     (3.11) 

% 

T 


fsin  W,t  -  sin  W0t] 


then  the  terms  S;k  (i,k=0,l)  in  Equations  3.8  and  3.9,   which 
are  a  func 
as  follows 


are  a  function  of  the  jammer  waveform  n.(t),  can  be  computed 


S„  = 


(S,,S),  +  (n.,S),|   + 

0 


(S,  ,C),  +  Cn.  ,c), 


a.  t 


and 


.2 


So,      [(S.,S)0  +  (n.,C)0j   +   [(S,,C)0  +  (n.,c)J 


P„,T 


(S0,S)0  +  (nk,C)c 


cos  4>  +  JiX 


(Sn ,C)  +  (n. ,C)fl 

•J 


4- 


f(S0,S),  +  (n  S)J   +   [CS0,C)(  +  (n.,C)J 


4- 
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Thus  the  probability  of  error  of  the  unmodified  receiver 
(Equation  3.10)  can  be  expressed  in  terms  of  SNR  and  JSR 
only  (defined  by  Equations  2.10  and  2.11)  as  follows 


V  f  x  -n 


^jHQ  (/f,/f>  -q  c/p^jd. 


(3.12) 


1 


Q  ( 


2  ,\T7 


0  -  Q 


/¥>}** 


where 


q(ix-     SNR  (2  +  2  J  2 JSR  cos  9  +JSP) 
oLo-     SNR  (2-2  7  2JSR  cos  <|>  +  JSR) 

$m=  p.o  =  snr  JSR 


Receiver  performance  can   now  be  evaluated  as   a  function  of 
SNR  for  fixed  values  of  JSR. 

To  provide  further  insight  into  the  performance  of  the 
modified  incoherent  receiver,  the  effect  of  varying  q/^  can 
be  analyzed  via  computation  of  the  probability  of  error  Pe 
from  Equations  3.8  and  3.9.  In  terms  of  SNR  and  JSR,  Pe 
becomes 


=  PCH.)   Pe,    +  P(H0)   ?Qo 

=  p(h.)  c,[ i  -  .gL-  q  (/cz^r^r,)  dej 

+  P(H|)   CAt7l\^°^  ^P^V^O  d0] 

•    +  p(hq)  c^i  ~\y  <Jc7Zo ,Jcfr0)  dtp 


+  P(H„)   C 


12 


^  j       Q   (£^0    JCiotoo)   d<p 


(3.13) 


42 


where 


C  , 


(!-*)" 


C 


The  effect  of  a   jammer  that  uses  a  single  tone   only  to  jam 
either  'mark'  or  'space'  signals,  can  be  analyzed  by  setting 


njCt)  =/V#-'inH't 


or 


n-(t)  =/pt17-4-£ 


T 


(3.14) 


If  the  first  choice  is  used,  we  obtain 


(3.15) 


?e    =  P(H,)  C, 
+  P(H.)  C- 


1  - 


J 

ZTl*o 


Z7L 


I  Q  (Jc,  o(u   Ju$d\  )  da] 
£q  (Jcx &',    Jc,  ad,'  )  de] 
+  P(H0)  C2  [  1  -  Jxffi  (JcTTTo    ,/iTpT  )  d<|>] 

[  ^/*q  (v^?r  ,y^ZT)  d«p] 


+  P(H0)   C 


where 


o/||=   2  SNP.  (1  +   2/TSR  cos  9  +  JSR) 


#= 

0 

</»» = 

2 

SMR 

(3p  = 

2 

SNR 

JSR 
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and  if  the  second  choice  is  used,  we  obtain 


Pe  =  P(H,)  C 


[  x  "  -^rlQ  {^Co  Jc*fi*  )  de 

P(H,)  C2       -±^  fcfwcip,',    JCit,.   )  d6 

1  -  -j^jQ  (Jc^]    JCifio',    )  ^9 


+  P(HJ  C 


+  P(He)  C, 


gTfjQ  (/^.^  ,/c777)  dcp_ 


The  effect  of  varying  the  threshold  oC  on  receiver  perform- 
ance for  the  jammer  waveform  of  Equation  3 . 14  is  discussed 
in  Chapter  5. 


C.   FREQUENCY  MODULATION  SWEEP  JAMMING 

In  certain  situations,  the  need  to  jam  a  certain 
frequency  band  rather  than  a  discrete  set  of  frequencies 
using  tone  jammers  may  arise.  Therefore  an  FM  sweep  jammer 
will  be  proposed,  analyzed,  and  its  effect  on  the  incoherent 
FSK  receiver  investigated  as  a  function  of  the  number  of 
times  the  jammer  sweeps  the  signal  band  during  the  signaling 
interval  [0,T].  The  mathematical  model  used  for  an  FM  sweep 
jammer  is 


n.(t)  = 
J 


\fi-*i*  [ 


^t  +  K^  fa  cos  W.t  dt]        0^t£T 
After    earring    out    the    integration   we    obtain 


n.(t)   = 
j 


%Phsini 


Wet  +   A  sin  W.t  + 


£ 
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where   p  =  K_ptf./Wj  and   g  is  a  deterministic   phase  angl 
The  instantaneous  jammer  frequency  is 


e . 


Wx(t)  =  Ws  +  o   W  cos  w#t 
r  J      J 
so  that  the   instantaneous  jammer  frequency  covers   the  band 

(Ws  -  ^W-  ,   Ws  +  £  Wj  )  .    Assume  that  W5T  =  2  7C 1  and  W.  T  = 

27Lk,   where  1   and  k  are  integers.    In   order  to  determine 

receiver  performance  in  the  presence   of  such  a  jammer,   the 

parameters  S;K  (i,k   =  0,1)   which  are  a   function  of  n.(t), 

o 

must  be  determined.   Before  so  doing,  we  evaluate 

(nJ>S)*   =  JjX^J^-sin   (  Wst  +  ^sin  W.t)  sin  Wkt  dt     (3.16) 


K  =  0  ,  1 

where  since  the  fixed  phase  &  represents  a  time  delay  only, 
it  is  set  to  zero  for  computational  ease.  Equation  3.16  can 
be  rewritten  as  follows 


(n.,S)K 
j 


T 


COS 


(Wc,  -  wK)t  +  P  sin  VLt 


-  cos  f(Ws+  WK)t  +  p  sin  W.t_ 


dt 


(3.17) 


By  using  the  well-known  Bessel  function  coefficient  expan- 
sion for  each  cosine  term  in  Equation  3.17,  the  following  is 
obtained 


(n.,s)K 


Z.1 


stJO^+wQ-hiw.It 


Si^^-^-f-rnalT       (3.18) 
(u^-u/kj  +  niu. 

T  — —  *J 


^s  +  w  j  t  n  vu 
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Since  the  FSK  signal  covers  approximately  the  band  (Wo  - 
4/t/T  ,  W,  +  47C/T)  (assuming  Wt  >  W0),  it  is  reasonable  to 
set 


Ws=j(W,  +  w0) 


(3.19) 


so  that  the  instantaneous  jammer  frequency  band  (W5  -  Q  W  .  , 
\i      +  Q  W«  )  completely  covers  the  signal  band  provided  that 

4-7C 


ws  -  £W.  =  w0  -  -£t£  j    ws  +  (5W.   =  W 


This  means  that 


p.  =-!-(  w,  -  w„)  +  4£ 


must  be  satisfied.   Since  for   FSK  signaling  we  have  assumed 
that 


(  W,  -  w0)  T  =  17C 


(W,  +  wD)  T 


m7TL 


where  1  and  m  are  integers,   and  W*T  =  27Ck,   it  is  apparent 
that 


(8=  -gfc-[i.  +  1/2] 


Note  now  that  the  integer  k  determines  the  number  of  times 
the  jamming  waveform  will  sweep  the  signal  band  in  one  bit 
interval.  Thus  from  Equation  3.18  and  Equation  3.19,  we 
have 


(n.,S), 

J 


r  srnf^Cu;.-u;()+tiu;.jT 

2"(3^  +  U)+  7)U/.   J 


which  becomes 
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p 


(n.,S), 


I  T 


Z. 


Z^(  (3 


Tl-.-oO 


Sin  7c(^K-Mn-fi/£) 


(3.20) 


and 


(n.,S)0 
J 


^_  £  Jfl(  p  )("    Sin[T(^-^)+WjT 


which  becomes 


(n.,S). 
j 


Bui 


g. 


£j,(p 


-«  =  -■> 


(3.21) 


If  1  is  an  even  integer,  then  Equation  3.20  and  Equation 
3.21  at  most  contain  two  terms  respectively.  These  terms  can 
exist  only  if  the  argument  of  all  SINC  functions  is  zero. 
Therefore  Equation  3.20  and  Equation  3.21  can  be  simplified 
to  yield 

(npS),  -~l^r-  [Jt,,(  P   -  J*a  <  (3)] 
where  n,  =  l/4k,  n^  =  -  (m  +  1/2)  /  2k  and 

-  (m  -  1/2)  /  2k. 


(n.,S)( 
j 


where  n: 


a 

l/4k,  n^ 


In  view  of  practical  communication  system  constraints, 
the  integer  m  will  typically  be  much  larger  than  1  because  m 
and  1  represent  the  sum  and  difference  of  the  signaling 
frequencies   respectively.     In   particular    we   see   that 
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JK(^)<<   1  for   ^  <<      k.    Hence  the   terms  J^  (  B   )    and 
J-n   (  6>  )  can  be  neglected  so  that 


(n.,S), 


JW   J-n,  <  (»  ) 


(3.22) 


and 


Cn.,S)0  =/M  Jl|j(  p  ) 


(3.23) 


The  other  parameters  to  be  considered  are 

(n.,C)K  =  I  Jl^.p=rsin   (W6t  +  (3  sin  W.t)  cos  WKt  dt 

where  0  =  Kr A  /  W- .  Using  again  the  Bessel  function  coef- 
ficient expansion  and  appropriate  trigonometric  identities, 
we  obtain 


(n.,C) 


-f 


-f(2nK-Je/2) 


(3.24) 


and 


CO 


(n.  ,  C)0  = 

J 


"  sin2'-g-(2nk  +  1/2)      (3.25) 


ti  =  -cp 


(2nk  +  1/2) 


sin*-^(2nk  +  m  -1/2) 
(2nk  +  m  -1/2) 
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We  note  that  all  the  terms  in  Equation  3 . 24  and  Equation 
3.25  take  on  the  form  of  a  SINC  function  times  a  sine  func- 
tion.  We  can  readily  show  that  for  even  values  of  m  and  1, 


Cn.,C),  =  Cn.,C)fl=  o 


(3.26) 


Using  the  mathematical  forms  given  by  Equations  3.22,  3.23 
and  3.26,  the  performance  of  the  receiver  with  FM  sweep 
jamming  is  obtained  via  the  use  of  Equations  3.8  and  3.9  and 
evaluated  from  the  expression 


p4  =pch.)  c.ti.^JQcy^  Jyzr^)de] 

+  pch.)  c2  [-^JTc/cT^  JZZ  )  de] 

+  p(Ho)  Ci  [-he  ${J^~°  ./gZT)  V 


(3.27) 


where 


dr  2  SNR  [  1  +  2/jSRj^  (  £  )  cos9  +  JSR  j-  (  p  }] 
oL=  2  SNR  [  1  +  27jSR  J^C  p  )  coscp  +  JSR  j^(  ^  j 
(&,=    (3,.  =   2  SNR  JSR  J^(fi) 
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C4  and  Cg,  have  been  defined  in  Equation  3.13.  As  before,  in 
this  equation  SNR  and  JSR  are  defined  by  Equations  2.10  and 
2.11  respectively. 

The  performance  results   to  be  presented  are  a  function 

of  SNR,   JSR,   and  the  number  of  times  the  jammer  sweeps  the 

signal  band.  This  is  discussed  in  greater  detail  in  Chapter 
5. 
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IV.  ANALYSIS  OF  FSK  IN  THE  PRESENCE  OF  JAMMING  AND  FADING 

A.   GENERAL 

In  certain  propagation  media,  the  received  signal  trans- 
mitted via  a  free  space  channel  is  often  subject  to  fading, 
a  phenomenon  caused  by  multipath  propagation  and  the  equiva- 
lent addition  of  random  phasors.  The  vector  sum  signal  will 
have  an  envelope  which  changes  with  time  resulting  in  an 
effect  known  as  fading. 

The  fading  signal  model  often  utilized  is  that  of  the 
nonselective,  slow  fading,  Rayleigh  distributed  signal 
amplitude  where  it  is  assumed  that  the  amplitude,  while 
random,  remains  constant  over  the  time  interval  [0,T].  Thus 
the  received  envelope  is  now  random  with  probability  density 
function 

fA(a)  =  _±£-exp(- Jl.)        a>  0 
A         E        E 

where  £  =  E{A2}  denotes  the  mean  squared  value  of  the  signal 
amplitude . 

The  optimum  (minimum  probability  of  error)  receiver 
structure  in  the  presence  of  additive  white  Gaussian  noise 
(however  with  no  jamming  present)  is  the  same  as  that  for 
noncoherent  nonfading  FSK  and  its  structure  is  therefore 
given  by  the  receiver  of  Figure  5. 

B.   ANALYSIS  WITH  NEAR  OPTIMUM  JAMMER 

We  shall  analyze  the  effect  of  Rayleigh  fading  on  an 
incoherent  FSK  receiver  operating  in  the  presence  of  a 
jammer  in  addition  to  the  additive  white  Gaussian  noise. 
Thus,  the  signals  at  the  front  end  of  the  receiver  are 
either 
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H,    :     r(t)    =  A  sin  (W,  t  +0)    +  n(t)   +  n.(t)        0  ^t  ^T 

J 


or 


H 


>(t)   =  B  sin  (VJ0t  +Cf>)   +  n(t)   +  n.(t)     O^t^T 


where  the   amplitudes  A   and  B   are  independent   identically 
distributed  random  variables  having  density  functions 


fA(a)  =  -^bexpC- 
fB(b)  =  -fcexpC--^) 


2  A, 


_b_ 


a  >0 


b>:0 


with  E{A}  =  E{B}  =  / -g  A0  and  E{A2}  =  E(B2}  = 


2A; 


Note  that  FSK  with  fading  is  the  same  as  what  was  previ- 
ously analyzed  (FSK  without  fading)  except  that  now  the 
signal  amplitudes  are  random  varibles.  Therefore  for  a 
fixed  value  of  A,  the  error  probability  when  a  "mark"  signal 
has  been  transmitted  is  the  same  as  that  of  the  noncoherent 
(nonfading)  FSK  case.  That  is,  we  can  use  the  results  of 
the  preceding  chapter  and  in  particular  make  use  of  Equation 
3.7,  to  obtain  (for  the  special  case  of  ^  =  1/2) 

2 


P. (A)  = 


2i\ 


(&■ 


Q(  _^£L,x)  x  exp 


,   xl+(^) 
xp(-     \  a  / 


) 


fS„ 


I.  (x 


')  dx 


de 


where  the  dependence  on   A  is  due  to  the  fact   that  the  term 
S 1 1   defined  previously  depends  on  A.   Then,  Pg   becomes 


P    = 


pe(A)    P(A)    dA 


1 
2tt 


faCi 


Q(     ,x)x    exp(- ) 


(4.1) 


■>0 
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I*  (x 


S.. 


)  dx   d 


4sexp(-  T^l)  dA 

i_i        ^         -  IX 


£.At 


o        L   Jo 


Q(^S,x)  x  exp(--|) 


~K  exp 


(  ~°^)-A 


2  A 


Io(xZp  dA  i  dx]  d8 


where  the  last  equality  has  been  obtained  by  interchanging 
the  order  of  integration.  This  equation  does  not  appear  to 
be  readily  simplif iable .     If  we  take  a  closer   look  at  the 


innermost  integral  in  Equation  4.1, 
evaluated  using 


I    (a    t)    exp(-p2t2)    dt 


that    integral   might    be 


v 


v 


1 IS. exp( Tl    F    fJi=H 

2pur(vU)  P      *p"lFl(    2 


-a 


+1    ,v+1»    Up2) 


where  Iv  (•)  is  the  modified  Bessel  function  of  the  first 
kind  of  order  v,  (•)  is  the  Gamma  function,  and 
,F,  (X,,X2;X3)  is  the  Confluent  Hypergeometric  function.  In 
the  special  case  of  X  ,  =0,  the  .  Confluent  Hypergoemetric 
function  ,F,  (O.X^jX^)  =  1  so  that  the  integral  in  question 
can  be  simplified  in  this  case.  However,  since  the  term  Sn 
includes  the  integration  factor  A,  we  may  be  able  to  calcu- 
late and  express  Equation  4 . 1  in  a  simpler  form  using  the 
formula  above  for  the  case  in   which  the  variable  A  could  be 


separated  out  in  the  term  /Sn   for  a  given  jammer  waveform. 

The  same  arguments  apply  to  obtaining  the  probability  of 
error  under  the  hypotheses  H0 .  Thus  Pe  can  be  shown  to  be 
given  by 

(4.2) 


'co     . —  2 

q(  l|i°,x)x  exp  (--j-  ) 


^° 


Az  exp 


O 
A  , z   ? 

(S^fs7o)-AZ 

.   Q  4 


I0  (x  -^2-)  dA   dx 
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d* 


Using  Equations  4.1  and  4.2,  we  can  evaluate  the  performance 
of  receiver  in  a  Rayleigh  fading  environment  from 


Pe  =  -T(Pel  +  Peo) 


assuming  that  the  two  signals  are  equally  likely  to  be  sent. 
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V.  DESCRIPTION  OF  GRAPHICAL  RESULTS 

A.  GENERAL 

In  this  chapter,  the  analytical  results  of  the  previous 
chapters  are  now  presented  via  graphical  means  based  on  the 
derived  mathematical  expressions  for  receiver  probability  of 
error . 

The  plots  presented  display  the  receiver  probability  of 
error  (Pe)  as  a  function  of  SNR  for  the  various  jammer  wave- 
forms previously  considered  for  a  set  of  JSR  values. 

In  each  plot,  the  case  JSR  =  0  has  been  included  in 
order  to  provide  the  basis  for  comparisons  of  the  jammer 
effectiveness  on  the  receiver  performamce  as  it  relates  to 
additive  white  Gaussian  noise  only  interference. 

B.  ASK  (ON  -  OFF  KEYING) 

The  graphical  results  for  the  incoherent  ASK  receiver 
performamce  are  presented  first.  These  plots  correspond  to 
numerical  evaluations  of  Equations  2.5  and  2.7  for  equally 
likely  hypothesis,   that  is,   P(H,)   and  P(H0)   are  equal  to 

1/2. 

The  plot  of  Pe  for  ASK  modulation  is  shown  in  Figure  6.7 
as  a  function  of  SNR  for  fixed  values  of  JSR,  using  a  jammer 
as  specified  in  Equation  2.6  .  Figure  6.7  clearly  shows  the 
'break  point'  phenomenon  in  which  if  JSR  increases  beyond  a 
certain  value  (0.25  in  this  figure),  Pe  increases  with 
increasing  SNR.  From  this  figure,  one  can  observe  that  16.0 
db  of  SNR  is  required  to  obtain  a  Pe  of  io~5  without  jammer, 
i.e. ,  at  a  JSR  value  of  0.0.  In  comparison,  it  takes  23.5 
db  of  SNR  to  obtain  the  same  Pe  for  a  JSR  value  of  0.1. 
Thus,  in  the  presence  of  a  jammer  with  a  JSR  value  below  the 
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break  point,  we  need  a  larger  SNR  in  order  to  obtain  the 
same  performance  level  of  a  receiver  operating  without  a 
jammer  interference.  However,  in  a  jamming  environment,  a 
JSR  value  above  the  break  point  produces  a  Pe  which 
increases  with  increasing  SNR.  In  fact,  there  is  no  value 
of  SNR  that  can  produce  Pe  of  10~;  for  a  JSR  above  the  break 
point.  For  the  case  of  ASK  modulation,  the  break  point 
occurs  at  a  value  of  JSR  which  is  approximately  0.25,  as 
obtained  from  Equations  2.8  and  2.9. 

Figure  6.8  shows  the  variable  thresholding  effect  on  the 
jamming  situation  with  JSR  =  0.3  beyond  the  break  point  (JSR 
=  0.25).  Instead  of  using  the  fixed  threshold  as  given  by 
Equation  2.5,  the  variation  of  the  receiver  threshold 
obtained  from  Equation  2.5  by  changing  the  value  of  R  ,  can 
reduce  the  jamming  effect  over  a  restricted  range  of  SNR 
values  as  shown.  However,  since  as  shown  in  Equation  2.13 
the  variation  of  the  value  R  does  not  significantly  affect 
the  threshold  value,  the  variation  over  a  wide  range  of 
values  of  R  does  not  result  in  a  significant  change  in  Pe . 

C.   FSK  WITH  TONE  JAMMER 

This  section  presents  graphical  results  pertaining  to 
jamming  effects  on  FSK  modulation  with  a  single  tone  jammer 
acting  against  one  of  the  two  channels  and  a  jammer 
consisting  of  two  different  tones  acting  against  both  chan- 
nels simultaneously. 

Figure  6.9  corresponds  to  the  performance  of  the  optimum 
FSK  receiver  in  which  pt  -  1/2.  Equation  3.10  or  Equation 
3.12  is  used  to  evaluate  performance  with  the  near  optimum 
jammer  specified  in  Equation  3.11  .  This  jammer  waveform 
can  be  thought  of  as  'mark'  and  'space'  channel  jamming. 
Figure  6.9  shows  a  similar  result  to  that  found  in  the  ASK 
case  except  that  the  breakpoint  occurs   at  a  higher  value  of 
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JSR  than  that   found  for  ASK.    This  breakpoint   occurs  at  a 

JSR  somewhere   between  0.5   and  1.0  as   shown  in   this  plot. 

From  this  figure  it  can  be  noted  that  13.5  db  of  SNR  is  need 

to  obtain  a  Pe  of  ]_0~5  f°r  a   JSR  value  of  0.0,   but  the  same 

Pe  is  obtained  by  increasing  the  SNR  to  16.5  db  for  a  JSR  of 

0.1.     This  demonstrates   that   relatively   low  JSR   values 

require   a  significant   SNR   boost  in   order   to  maintain   a 

certain  desired  Pe  value. 

As  shown  in  Figures  6.7  and   6.9,   comparison  of  ASK  and 

FSK  modulation  reveals  that  FSK   is  somewhat  less  vulnerable 

to  jamming.   However   it  must  be  remembered   that  the  jammer 

waveform  n.(t)  used  in  each  case  is  different. 
j. 

The  effect  of  the  single  tone  jammer  on  the  optimum  FSK 
receiver  is  presented  in  Figure  6.10  which  corresponds  to 
the  evaluation  of  Equation  3.15  .  The  single  channel 
jamming  on  either  the  'mark'  or  the  'space'  channel  has  the 
same  effect  insofar  as  single  tone  jamming  is  concerned. 
Therefore  the  effect  of  'mark'  channel  jamming  only  is  eval- 
uated and  plotted.  Note  that  in  Figure  6.10,  19.5  db  of  SNR 
is  required  to  obtain  a  Pe  of  iq-5  as  compared  to  an  SNR  of 
24.5  db  in  Figure  6.9  for  a  JSR  value  of  0.3  with  Pe  =  iq~5 
also . 

As  expected,  Figures  9  and  10  demonstrate  the  fact  that 
single  channel  jamming  is  less  effective  than  simultaneous 
jamming  of  'space'  and  'mark'  channels  with  a  near  optimum 
j  ammer . 

The  effect  of  a  variable  threshold  on  FSK  will  be 
considered  by  changing  the  value  of  q£  in  the  modified 
receiver  shown  in  Figure  6.6  .  For  a  value  of  qC  other  than 
1/2,  the  simultaneous  jamming  of  'mark'  and  'space'  channel 
results  in  a  compensation  of  the  other  channel  such  that  the 
jamming  effect  remains  the  same  as  in  the  case  of  qC  =  1/2. 
In  other  words,  it  is  difficult  to  reduce  the  near  optimum 
jamming  effect  by  means  of  a  varying  the  threshold. 
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On  the  other  hand,  when  one  channel  jamming  is  applied 
to  the  receiver,  the  jamming  effect  can  be  reduced  by 
adjusting  the  threshold  with  increasing  SNR  as  shown  in 
Figure  6.11,  part icularily  for  JSR  of  1.0  and  various  values 
of  #( .  A  moment's  reflection  will  reveal  that  the  'mark' 
channel  jammer  increases  the  output  power  level  of  the  upper 
envelope  detector  (see  Figure  6.6)  so  that  when  the  'space' 
signal  is  sent,  the  error  increases.  Therefore,  for  this 
type  of  error,  under  the  assumption  that  the  'space'  signal 
has  been  transmitted,  Pe  can  be  reduced  by  lowering  the 
level  of  the  output  of  the  multiplier  by  using  an  appro- 
priate value  of  £>( . 

D.   FSK  WITH  FM  JAMMER 

This  section  presents  the  effect  of  an  FM  sweep  jammer 
using  sinusodal  modulation  on  noncoherent  FSK  signaling.  The 
FM  jammer  was  designed  to  sweep  the  bandwidth  occupied  by 
the  signal  several  times  during  a  bit  interval.  Thus  the 
jammer  effectiveness  was  investigated  as  a  function  of  the 
number  of  times  the  jammer  sweeps  over  the  band  of  the 
signal  during  one  bit-time  interval. 

Figure  6.12  shows  the  result  for  one  sweep  of  the  jammer 
per  bit  interval.  Figure  6.13  shows  the  result  of 
increasing  the  sweeping  to  two  sweeps  per  bit  interval. 
These  plots  show  that  the  FM  sweep  jammer  can  be  more  effec- 
tive by  increasing  the  number  of  times  of  jammer  sweeping 
during  a  bit  interval.  This  can  be  expected  from  the  results 
obtained  in  Equation  3.27  .  In  Figure  6.12,  in  order  to 
obtain  a  Pe  of  io~5  ,  16.0  db  of  SNR  is  required  for  JSR  value 
of  0.3,  but  in  Figure  6.13  which  corresponds  to  twice  the 
jammer  sweeping,  the  same  Pe  can  be  obtained  by  increasing 
the  SNR  value  to  19.5  db  for  the  same  JSR  value  of  0.3.  In 
comparison  to   the  previous  case   of  FSK  with   tone  jamming, 


for  the  same  jamming  environment  (i.e.,  JSR  value  of  0.3) 
Figures  6.9  and  6.10  show  that  24.5  db  and  19.5  db  of  SNR 
are  required   respectively  in   order  to  get   the  same   Pe  of 

10"5. 

These  different  requirements  of  SNR  value  for  various 
jammer  waveforms  show  that  the  FM  sweep  jammer  can  be  effec- 
tive but  in  general  is  not  as  effective  as  the  near  optimum 
jammer.  Note  that  from  a  practical  point  of  view,  the  added 
complexity  of  FM  jammer  waveform  may  make  it  an  unlikely 
candidate  for  a  replacement  of  the  near  optimum  jammer. 
However,  one  advantage  the  FM  sweep  jammer  has  over  the  near 
optimum  jammer  is  that  the  former  can  spread  its  power  over 
a  large  bandwidth  easily  and  therefore  is  more  effective 
than  the  latter  in  the  case  of  lack  of  exact  information  or 
knowledge  about  the  signal  carrier  frequency. 
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VI.  CONCLUSIONS 

The  familiar  model  in  which  Gaussian  noise  is  the  total 
interference  is  not  adequate  when  jamming  or  interference 
signals  are  present  in  the  transmission  environment.  This 
thesis  has  analyzed  the  effect  of  various  deterministic 
jammer  waveforms  in  terms  of  probability  of  error  on  binary 
incoherent  receivers  operating  in  the  presence  of  noise. 

From  the  jammer  point  of  view,  the  goal  is  to  cause  the 
maximum  possible  error  to  the  various  receivers  while  making 
efficient  use  of  its  available  power  (i.e.,  with  fixed 
jammer  power).  For  coherent  receivers,  it  was  proved  that 
the  optimum  jammer  waveform  is  made  of  a  deterministic 
signal  proportional  to  the  difference  of  the  binary  signals 
used  to  carry  the  digital  information.  This  thesis  has 
demonstrated  that  those  optimum  jammers  derived  for  coherent 
receivers  perform  their  function  as  near  optimum  jammers 
satisfactorily  against  incoherent  receivers. 

Therefore,  these  nea-r  optimum  jammers  can  be  concluded 
to  be  one  of  the  most  attractive  candidates  for  efficient 
jamming  of  binary  incoherent  communication  systems.  An 
optimum  jammer  has  not  been  derived  or  analyzed  because  the 
complexity  of  the  expression  for  Pe  makes  it  very  difficult 
if  not  impossible  to  derive  the  optimum  jammer  waveform  in 
closed  form.  Other  jamming  waveforms  such  as  single  channel 
jamming  and  FM  sweep  jamming  showed  its  inferiority  in 
comparison  to  near  optimum  jammer  waveforms  and  their  effi- 
ciency can  be  reduced  by  means  of  appropriate  variations  of 
the  receiver  threshold. 
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Figure  6.1    Correlation  Receiver  for  Binary  Signal 
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Figure    6.2         Quadrature    Receiver    for   ASK 
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Figure'  6.3         Alternate    Form    of    Quadrature    Recei 
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Figure  6.4    Incoherent  Matched  Filter  Receiver. 
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Figure  6.5'   Incoherent  Frequency  Shift  Keying  (FSK)  Receiver, 
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Figure  6.6    Modified  Incoherent  FSK  Receiver 
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APPENDIX  A 
DIGITAL  COMPUTER  IMPLEMENTATION  OF  THE  MAR CUM  Q-  FUNCTION 

The  Marcum  Q- function  occurs  frequently  in  communication 
problems  involving  incoherent  detection  of  signals  with 
single  or  multiple  observations  in  the  presence  of  noise  and 
jamming.  So  it  is  often  necessary  to  compute  values  for  the 
Marcum  Q-function  which  is  defined  by 

Q(  a  ,  6  )  =    V  exp(-  ?—±-    )    IQ(o  V)  dV 
^  2 

=  1  -    V  exp(-   V  +a   )  IQ(a  V)  dV 


where  1^  ( •  )  is  the  modified  Bessel  function  of  zero  order. 
It  is  noted  that  the  integrand  f  (  V  ,  pO  of  Equation  A.l  is 
the  Rician  density  function  which  is  sometimes  called  the 
generalized  Rayleigh  density  function.  The  normalized 
Rician  distribution  for  (y  2  =  1  is  shown  in  Figure  A.l  . 
For  large  values  of  q(  ,  the  Rician  density  function  f (V) 
can  be  approximated  by  the  normalized  Gaussian  density  func- 
tion with  mean  value  of  approximately  q^  .  This  can  be  justi- 
fied as  follows 

2   2 
f(V)  =  V  exp(--!L+2L  )  iQ(a  V)  dV 

v24.~2       aV 
~  V  expC-^-  )-?A-r 

f  ^2 

2  ( v-ct ) 


y2anV 


7iV      L    2    J  ^ 
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0  6r 


2   NO(S£  POwtR 


Figure  A.l    Rician  Density  Function 

_  X 


/27tx' 


and  ^V  ■  V1 


by  use  of  the  approximations  IQ (X) 
for  large  ^ . 

In  order  to  implement  the  computations  of  the  Marcum 
Q-function  using  a  digital  computer,  particular  interest 
must  be  focused  on  values  of  the  function  on  the  tails  of 
the  Rician  density  defined  by  the  integrand  in  Equation  A.l 
Specifically,  small  values  around  the  tails  of  the 
density  must  be  monitored  because  of  the  limitations  of  the 
digital  computer.  By  a  change  of  variables  (X  =  g^  V)  the 
Marcum  Q-function  defined  above  can  be  expressed  as 


Q(  a  ,  3  )  =  1  - 


x  exp 


-x 


where  f(X)  -Q  I0(X).  Since  the  digital  computer  can  treat 
the  values  in  the  range  of  y  such  that  -180.218  £  y  5 
174.643  for  the  function  eY  (  for  the  IBM  3033  ),  for 
digital  computer  implementation  of  the  Marcum  Q-function  the 
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integral  limits  need  to  be  adjusted  to  meet  the  conditions 
for  these  acceptable  values  in  the  computer.  Negligible 
areas  around  the  tails  of  the  density  beyond  certain  bounds 
must  be  discarded  without  significantly  increasing  overall 
error.  Therefore  from  Equation  A. 2,  the  integral  limits  can 
be  substituted  with  appropriate  bounds  such  that  exponent 
satisfies  the  condition  given  by 
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In  other  words,  the  variable  X  must  be   located  within  the 

range  p<  2  -  qL  /  3  60  <  X  <  ^  2  *  ^  y  3 6 0  -  For  a  value  of  X 
beyond  that  range  the  exponential  of  Equation  A. 2  is  too 
small  (less  than  Q  )   compared   to  other  computed  values 

so  that  the  area  outside  of  the  new  limits  can  be  neglected. 
Let   us  denote   UL   and  LL   the   upper   and  lower   limit 
respectively  defined  by 


UL  =  a   +  a/360" 


0  r-^— 

LL  =  a   -  a/360 


Then  we  can  consider  the  various  situations  case  by  case. 
If  the  lower  limit  LL  is  positive  for  certain  values  of^(, 
then  in  the  case  of  q(Q  <  LL ,  that  is,  &  <  ot  -  /3  60,  the 
integral  value  can  be  ignored  so  that  the  value  of  Q(^,p  ) 
can  be  assigned  to  be  one.  In  the  case  of  LL  <  p^B  £  UL, 
that  is,  qC  "  /360  <  Q  <  pt  +  /360,  the  Marcum  Q-function 
can  be  approximated  by 


HP 

Q(  a  ,  S  )  =  1  -  —i\    x  exp 


'LL 


,  x     .2 
(  —  -  a ) 


f(x)  dx 


In  the  case   of  M    >  UL,   that   is,   /$  >  oC   ♦  J   3  60,   the 
Marcum  Q- function  can  be  computed  from 


Q(  a  ,  6  )  =  1  -  — 


•UL 

x  exp 


'LL 


(  —  -  a) 


f(x)  dx 
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On  the  other  hand  if  the  lower  limit  LL  is  negative  for  some 
value  of  q£  ,  then  there  are  two  cases  to  be  considered.  In 
the  case  of  ^  <  UL,  that  is,  &  <  o(  +  /  360,  the  Marcum 
Q- function  can  be  substituted  with  th^  value  given  by 


46 


QC  a    ,  B  )  =  1  -  — 


x  exp 


a 


Finally  in  the  case  of  oLB    >  UL,  that  is,   B  ^  (/.    + 
the  value  of  Marcum  Q-function  can  be  computed  from 


Q(  a  ,  3  )  =  1  - 


a 


/X     .2 
(—  -  a) 

a  


f(x)  dx 


Here  it  has  been  assumed  that  the  function  f(X)  defined  by 
e  IQ(X)  does  not  impose  the  limitation  of  computation  on 
the  digital  computer  and  that  the  library  functions  for  the 
computation  of  e  Ij(X)  and  its  integration  with  desired 
accuracy  are  available  to  the  user. 
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